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QJ . Abstract 
^ ■ 

Charm and bottom transport coefficients in a medium constituted of light mesons, such as is formed in the hadronic phase of Heavy 
Ion Collisions, are obtained within an effective field theory approach implementing heavy quark symmetry and chiral symmetry 
breaking. Heavy flavor propagates in the medium as D/B and D' IB' degrees of freedom, and unitarization of the lowest order 
^ ' heavy-light meson amplitudes is used in order to reach high temperatures. The latter accounts for dynamically generated resonances 

' in isospin 1 /2 channels, a feature that leads to a more efficient heavy flavor diffussion. We discuss the temperature and momentum 

dependence of the friction and diffusion coefficients in a transpoit approach up to temperatures of about T 150 MeV, and provide 
estimates of the charm/bottom relaxation lengths and momentum loss. Implications for heavy meson spectrum observables in 
Heavy Ion Collisions are discussed. 
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, 1. Introduction 



The features of matter formed in Heavy Ion Collisions (HICs) have been a subject of great interest in the last 
decades. In this scenario, heavy-flavored hadrons play an essential role since they carry heavy quarks produced in the 
early stage of the collisions, unlike pions and kaons which can be produced in the thermal medium at later stages. 

It is worth noticing that the momentum spectra of charmed and bottomed mesons extracted from HICs undergo 
modifications due to their interactions with the hadronic medium. In this sense, the diffusion of heavy mesons in an 
equilibrium hadronic gas should be taken into account in order to compute realistic spectra, for instance, in transport 
approaches. Unlike other lighter systems, relatively Uttle attention has been devoted to study heavy-meson dynamics 
in the hadronic phase of a HIC. In addition, there is a considerable dispersion in existing results on the matter, as for 
instance concerning the value of the charm relaxation length at a given temperature 1 1, 2, 3, 4, 5, 6]. 

In this talk we report on recent progress in the calculation of the transport coefficients of charmed and bottomed 
mesons in a hot gas consisting of pions, kaons and rj mesons. In order to do this we provide a realistic determination of 
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the relevant heavy-light meson scattering ampUtudes in an effective field theory approach, exploiting chiral perturba- 
tion theory at next-to-leading order (NLO) together with constraints from heavy quark symmetry, and implementing 
exact unitarization in order to reach higher temperatures. The latter is a crucial point in our analysis, leading to dy- 
namically generated states in agreement with our current knowledge of the D and B spectra with a minimal number of 
parameters, and therefore to a more realistic (and theoretically consistent) determination of the transport coefficients. 



2. Modelling heavy-meson interactions in a light meson gas 

We calculate heavy-meson transport coefficients in a Fokker-Plank equation approach. The dynamical input is 
encoded in the scattering amplitudes of (long-lived) charmed and bottomed mesons with the octet of light pseudoscalar 
mesons. We assume that the density of heavy mesons is negligible and thus ignore collisions among themselves. For 
a detailed discussion of the model for charm interactions within the pion gas we refer to our previous work [Q. 

The amplitude for scattering off' a charm (bottom) quark in the light meson gas, at next-to-leading order (NLO) in 
the chiral expansion and leading order (LO) in the heavy quark expansion, irrespective of whether the heavy quark is 
in a D{B) or a D{By meson state, is given by 

^ = - m) + yr'^^ + 7T'^3(P2 ■ P4) + yr'^^ ^^P^ ■ /'2)(P3 ■ P4) + (Pl ■ Pa){P2 ' Pi)\ , (1) 

where C,- (/ - 0, 3) are channel-dependent numerical coefficients in isospin basis. The three (ffavor-dependent) 
parameters /i, are the low-energy constants (LECs) from the NLO chiral Lagrangian, which we constrain with the 
available experimental information on the heavy-meson spectrum (a list of the relevant states is accounted for in 
Tableffli. 

The chiral perturbative expansion is typically bound to work properly only at very low energies, and cannot 
describe the appearance of resonances in a given scattering channel. In addition, at high energies, perturbative cross 
sections violate Froissart bounds imposed by the unitarity of the S -matrix. We solve these limitations by implementing 
exact unitarity in the scattering amplitudes by solving the Bethe-Salpeter equation with the perturbative amplitudes as 
dynamical kernel, namely 

T{s)^-V{s)l{\-V{s)G{s)), (2) 

where V{s) is the S -wave projection of the scattering amplitude in Eq. ([1), and G{s) stands for the two-meson resolvent 
function (loop integral), conveniently regularized. 

The free parameters of the model include a subtraction constant from dimensional regularization of G(s), plus 
the three LECs. The subtraction constant is fixed by the resonance position in the isospin 1/2 channels. Note that 
these resonances are already generated by the LO in the interaction [{s - u) term in Eq. ([T]i], which is completely 
determined by chiral symmetry breaking, and the freedom provided by the LECs can be used to improve on the LO 
results. They have to be determined for the charm and bottom sector, although some of the LECs scale with the heavy 
meson mass and thus not all of them are independent. Note that the bottom sector is much less known experimentally: 
Particularly, the 5 -wave resonances have not yet been observed. Given the fact that the two D\ states are very close 
in energy, we assume that the same holds for the bottom sector, provided that heavy quark symmetry works. Also 
note that the D{B)n and D{B)*n amplitud es are related in the heavy quark limit just by exchanging mo{B) by mD(B)' 
in Eq. (|2]i. For the bottom sector, we obtain as our best determination ma, = 5534 MeV (Fb,, - 210 MeV) and 
mfi, - 5587 MeV (Fb, = 250 MeV) in good agreement with III IS Hi (cf. 101 for the charm sector). Threshold cross 
sections (or scattering lengths) are also obtained in good agreement with previous determinations. 



3. Heavy-meson transport coefficients 

The momentum-space distribution of bottomed mesons must relax via the Fokker- Planck equation (we refer to 
for details and assumptions). Let us label the momenta of the elastic collision between a heavy meson M^*' and a light 
meson (p as M'**(p) + 0(q) — > M'**(p - k) -H (f>{q + k). The evolution of the momentum distribution of the heavy meson 
due to its interaction with the isotropic mesonic gas is controlled by the drag (F) and diffiision (Fo,Fi) coefficients, 
written as 

F(p^) ^ Jdk w(p, k) ^ , Toip^) ^Ijdk w(p, k) 



k^- 



(kip-y 



ri(/) 



4/ 



dk w(p, k) 



ikiP') 



l\2 



(3) 
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Spin, J" 


D (QM) 


D (exp.) 


(M, r) MeV 


this work 


B (QM) 


B (exp.) 


(M, T) MeV 


this work 


1/2, 0+ 


Do 


D*(2400) 


2318,267 


2300, 350 


Bo 


7 


7 


5534,210 


1/2, 1 + 


Di 


Di(2430) 


2427, 384 


2390, 400 


By 


7 


7 


5587, 250 


3/2, 1 + 


Di 


Di(2420)° 


2421, 27 




Bi 


Bi(5721) 


5723, ? 




3/2, 2+ 


D2 


£);(2460) 


2466, 49 




Bi 


(5747) 


5743, 23 





Table 1. List of heavy meson states. Left side: D meson. Right side: B mesons. The spin 1 /2 doublets are the ones dynamically generated from 
unitarization of the S-wave interaction. Experimental data from fioll . Lack of experimental evidence is labeled with a question mark. 



where w(p, k) is the colHsion rate for a heavy meson with initial (final) momentum p (p - k), 

(4) 

Here /^(q) is the bath's distribution function, Mm^ stands for the Lorentz invariant heavy-light meson scattering 
matrix element, is the Goldstone boson isospin degeneracy (e.g. g„ = 3 for the pion), and^ denotes possible spin 
degrees of freedom. Since spin average is trivial in view of heavy quark symmetry, one just has 2 \Mm<p{s, t,x)\^ = 
Y^^(2i+i) 2/(2/ + l)|r'|^, in terms of the isospin averaged unitarized amplitude. 

In Fig. [T] we show the temperature and momentum dependence of the B-meson drag and diffusion coefficients. 
We observe an increase of factor 6-8 in the range T = 100- 180 MeV, namely the drag in a heavy-ion collision 
is considerably strengthened in the hotter stages, with significant interaction between heavy mesons and the thermal 
medium (and accordingly at larger momentum transfers). From F we can estimate the relaxation length of bottom 
quarks in the hadronic medium. In the pion gas at T = 150 MeV we find Ab{T = 150 MeV, p = I GeV) - j ^ 100 fm 
for a typical momentum of 1 GeV, to be compared to Ad - 40 fm for charmed mesons evaluated in the same approach 
101 . Therefore bottomed mesons barely relax during the lifetime of the hadron gas, unlike charm mesons that, while not 
relaxing completely, may loose a great deal of memory of the initial state. In this sense, bottomed mesons constitute an 
optimal system to characterize the early stages of a relativistic heavy-ion collision, in spite of their (detected) number 
being less abundant than for charm. Also note that some of the observed c quarks in the final state proceed from b 
quark weak decays in the micrometer range, and therefore they interacted as b quarks with the (hadronic) medium. 

We observe that, in the static limit (p — > 0), the drag coefficient for bottomed mesons is about 3 times smaller than 
for the charm case, whereas the diffiision coefficients are very similar for the two systems, c.f. Fig. [1] This scaling 

is the correct one as expected from non-relativistic kinetic theory, F ^ |crP -yj^i^ ™d Fq, Fi ^ |crP y/m^T (cr is 
the total cross section and P is the pressure of the gas). Such result constitutes a consistency test of our calculation 
and is also observed for the transport coefficients of heavy quarks beyond the critical temperature iullil. We have 



found that implementing unitarity to obtain realistic amplitudes for heavy meson scattering off the light gas plays an 
important role and has a strong impact in the final size and temperature dependence of the transport coefficients, as 
well as on their scaling properties with the heavy quark mass. For instance, the use ofperturbative cross sections 
overestimates scattering at high energies and leads to an unrealistically large drag force 1131. Alternatively, scattering 
lengths have also been used as dynamical input |6(], which underestimates difussion due to the i-channel enhancement 
of the interaction in the resonant channels. The unitarized scheme provides the most realistic result, coherently with 
our current knowledge of the heavy meson spectrum and with a controlled high energy behavior 

We analyze the contributions from the different species in the S U (3) gas for a bottomed quark in Fig. |2](left). As 
expected, the most relevant contribution comes from the pion gas. At T = 150 MeV, pions provide almost 90% of the 
total, while the next-to-leading contribution is provided by kaons and (mostly) anti-kaons. The relaxation length of 
bottomed mesons travelling with 1 GeV momentum in the full meson gas is Ab(T - 150 MeV, p - 1 GeV) - j - 
80 fm, c.f. Fig. |2](center). Still, this thermal relaxation time is greater than the lifetime of the hadron gas, reinforcing 
the idea that heavy quarks are optimal carriers of information of the phase transition upon exiting the interacting 
region. Our estimate of the spacial diffusion coefficient for charmed and bottomed mesons. Fig. |2] (right), together 
with estimates in the quark-gluon plasma phase, indicates that the relaxation time of heavy quarks has a minimum 
around the crossover 
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Figure 1. Momentum-space drag (F) and diifusion (Fo, Fi) coefficients as a function of temperature for several bottom quark momenta in the pion 
gas. The charm coefficients are shown for comparison at p = 0.1 GeV 0]. 




Figure 2. Left: Contributions to the Fq coefficient in the SU(3) meson gas at p = 0. 1 GeV. Center: Bottom relaxation length in the hadronic phase 
at p = 0.1 GeV as a function of the temperature. Right: Charm and bottom spacial diifusion coefficients below (this work) and above ([Tll[l2,[i3|l) 
ffie crossover 



Finally, a momentum loss per unit length of about 50 MeV/fm (70 MeV/fm) is obtained in our approach for a 
charmed (bottomed) meson propagating in the gas with a typical momentum of 1 GeV with respect to the rest frame of 
the medium. The results discussed above allow us to estimate the modification in heavy quark spectra due to diffusion 
in the hadronic phase. The observed distribution of D mesons, being different from the one after hadronization, may 
lead to a different determination of the freeze-out temperature from fits to the Botzmann exponential shape at low 
Pt- Observables like the nuclear suppression factor or the elliptic flow of heavy flavored mesons can be sensitive, 
in particular, to the sizable momentum dependence of the transport coefficients found in this work, which could be 
implemented in hydrodynamical simulations ll4 l. 
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